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Abstract. In this work we introduce an energy function in order to study fi- 
nite scale free graphs generated with different models. The energy distribution 
has a fractal pattern and presents log periodic oscillations for high energies. 
This oscillations are related to a discrete scale invariance of certain graphs, 
that is, there are preferred scaling ratios suggesting a hierarchical distribution 
of node degrees. On the other hand, small energies correspond to graphs with 
evenly distributed degrees. 



1. Introduction 

In the last few years a huge amount of research on graphs has been achieved 
in physics and mathematics, mainly due to the increasing importance of complex 
networks like the Internet and the World Wide Web, and the role of social networks 
in the propagation of diseases, infections, rumors and news, in both real and virtual 
environments such as populations or web-based communities. 

Despite the different origins, methods and perspectives in those works, certain 
concepts appeared and pervaded the network literature, like small worlds, cluster- 
ing, centrality, scale free degree distribution, and preferential attachment, among 
others, pQ, 0, [HJ- They suggest that simple underlying mechanisms could describe 
and explain the network growth and their emerging properties. To this end, several 
networks based on real data were analyzed (see 0], [S] and the references therein), 
and many theoretical models of network growth were proposed (see Table 3 in 

Also, there are still many questions and problems which remain unanswered. 
Perhaps the main one is the origin of the power law degree distribution of nodes 
in real networks. This distribution seems to be observed in several real networks 
(although some criticism appeared in the last year, [Jj, [H]). Roughly speaking, 
the power law distribution says that the number of nodes in a graph with k links 
for k > fco has a decay proportional to k~ a , for certain a > 1. Several theoretical 
models were proposed to explain this phenomenon. Among them, let us mention 
preferential attachment of Barabasi and Albert [H], edge redirection of Krapivsky 
and Redner and attachment to edges of Dorogovtsev, Mendes and Samukhin 

m 

Although much attention was paid in the theoretical literature in the thermo- 
dynamical limit (when time and number of nodes tend to oo) in these models, not 
too much is known for graphs with a bound in their size being an interesting 
exception) . Indeed, many of the computer-based simulations tend to concentrate 
on large graphs, while small ones usually get ignored. Yet, real networks are fi- 
nite, and in many cases they are very small (like networks of metabolic reactions, 
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proteins interactions, digital electronic circuits, sexual interactions or food webs). 
Our intention in this paper is to concentrate on these small graphs. Specifically, we 
concentrated on graphs with a number of nodes ranging from 50 to 100,000. On 
the other hand, we were able to generate large amounts of graphs for each model, 
which reduced the noise and statistical fluctuations, and helped see some otherwise 
hidden properties. 

Another question, which motivates our research, was the discrimination among 
the different graphs generated from these models and the analysis of their proper- 
ties. We may consider this question as a basic one concerning the internal topology 
of the network, and the quality of their connections. Let us note that different net- 
works with the same node distribution respond in distinct ways to edges or nodes 
removal depending on the connectivity between hubs, that is, if the nodes with high 
number of links are connected or not between them. This question was considered 
by Newman in |12) . who coined the term assortativity in order to describe this 
phenomenon. 

In order to study this question, we introduce the linking energy of a graph G, 
defined as 

E 2 (G) = J2\d(x)-d(y)\ 2 , 

where d(x) is the degree of the node x, i.e., its number of links, and the sum runs 
over the pairs of connected nodes. 

Our original intention was to compute the expected linking energy for certain 
scale-free network generating models, which enabled us to classify the graphs gener- 
ated for each models as either assortative or not. One of the first observations, was 
that the underlying mechanism responsible of networks generation in each model 
could generate both assortative and non-assortative networks. Another interesting 
fact is the correlation of the assortativity and the linking energy E 2 , and the fractal 
structure of this correlation. These results will be published in a separated work. 

However, we also obtained an unexpected result when we computed E 2 for the 
models mentioned above, since the distribution of energies presents log-periodic 
oscillations in the tail. This oscillatory behavior is more apparent and becomes 
amplified within models which tend to produce a low number of highly connected 
hubs. The oscillatory behavior is not present for random or regular graphs, even 
when they are small. 

Surprisingly, this shed some light on the first problem, the node distribution. 
The presence of log-periodic oscillations is associated to fractal complex dimen- 
sions, which are due to discrete scale invariance instead of a scale free similarity 
(dSj, El)- Indeed, this suggests a stronger order in the networks other than that 
predicted by the scale free models, since fractal complex dimensions are related 
to (mathematical) self similar fractals such as the triadic Cantor set. Hence, we 
may classify the graphs generated with any of those models as continuous or lacu- 
nary scale free. A similar dichotomy is present for fractals generated with iterated 
functions systems |15| . 

We conjecture that the nodes distribution itself must also exhibit the oscillations. 
Further work is required to settle this question. However, several real networks that 
appeared in the literature seem to show such oscillations, like the World Trade Web 
[Tfi| . or web-based communities (see Fig. 1, 5 and 8 in 53)- Also, two theoretical 
models of networks without growth were recently proposed in and > where 
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oscillations are apparent. The existence of oscillations on the degree distribution 
was related to the existence of different hierarchies among the nodes of the graph, 
see j2J where theoretical models were proposed to study this kind of networks. 
Hence, our results show that classical models based on preferential attachment or 
edge redirection can be used to generate non homogeneous networks. 



2. Main Results 

2.1. The linking energy. Let G be any graph, and let us denote x ~ y whenever 
there exists a link between the nodes x and y. Let d(x) be the degree of x, that is, 
the number of links attached to x. 

We define the linking energy of G as 

E 2 (G) = J2\d(x)-d(y)\ 2 . 

This is a discrete version of the energy which arises for the tension of a system of 
strings attached at the nodes, with heights proportional to the number of connec- 
tions. For regular networks, it coincides with the variational form of the Dirichlet 
energy of a discrete Laplacian. 

Let us note that the energy spectra for graphs with N nodes is finite. For a tree 
on N nodes, the maximum value of the energy is attained for the central tree where 
all the nodes are connected to one of them, with an energy of 

(N- l)(N -2) 2 ~ TV 3 . 

Let us note as well that the energy is an even number: 

e 2 {g) = j2d x -dy = J2 dx+d y = J2 d ^ = J2 dx = J2 2 =° ( mod 2 ) 

Xr^y Xr^y X X X^y 

For regular lattices, such as the triangular, square or hexagonal lattice, this 
energy is proportional to the number of boundary nodes, since the degree is the 
same for all the internal nodes. Also, a random (Erdos-Renyi) graph has lower 
energy than a scale free one due to the small number of hubs in it. 

For networks with power law degree distribution, this degree sequence is not 
enough to characterize several topological and dynamical properties. Graphs with 
the same degree sequence behave differently, depending on whether or not the hubs 
are connected among them, a problem considered by Newman who introduced the 
assortativity to measure the connectivity between hubs. The linking energy could 
be used in order to study this problem. 

2.2. E% in a random (Erdos-Renyi) graph. Let G be a graph in the Erdos- 
Renyi model, with d nodes. Here, two nodes connect each other with probability 
p. Let us fix an order for the nodes, say x\, ■ ■ ■ ,Xd, and for 1 < i < j < d let 
Xij = #{fc ~ i, k 7^ j} and let Y.- L j = #{fc ~ j, fc ^ i}. In other words, JQj is the 
degree of the node i disregarding a possible connection to node j, and similarly for 
Yij. Then, and Yy are binomial, with probability 

piXij =k)= P {x n = k) = - 2 y (i - P ) d - 2 - k . 
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Thus, the mean and the variance are given by 

< X l3 >=p(d-2), 

a 2 Xi . =< X% >-< X {j > 2 = (d - 2)p(l - p). 
Furthermore, the variables Xij and Y{j are independent. Hence, 
< {Xij - Y iS f > = 2<Xf j >-2< X id > 2 
= 2(d-2)p(l-p). 
Now, the mean value of E2 is given by 

<E 2 >=Y,P<{X i j-Y ij f > 

'[jp-2(d-2)p(l-p) 

= d{d-l)(d-2)p 2 (l-p). 

For other models of graphs, the distribution of energies gives better information 
than the mean value of E2 . 

2.3. Models of graph generation. We explain now the studied models. All of 
them consist of the following steps: 

• One begins with a fixed seed of nodes and links. 

• At each step, a new node is added to the graph, and 

• a fixed number of links emanate from the new node to the existing ones. 
These models differ in the way the targets for the new nodes are chosen. We 

considered the following models: 

BA Preferential attachment (or linear P.A.): one target node is selected at 

random, with a probability proportional to the degree of the nodes. 
KR Edge redirection: one target node is selected at random with uniform 

probability. However, with probability 1 — r, this target is changed by the 

node it points to. 

DMS Attach to edges: one link is selected at random with uniform probability 
(among the links). The new node connects to both ends of the chosen link. 
The labels are put after the authors of these models, see [5], ^UJi [TTj . 

2.4. Numerical Experiments. We computed the energy E% for several graphs 
constructed according to these models. We present in this subsection some of these 
computations. For the KR model, we used r = 0.5, which gives a power law with 
exponent 3 for the degree distribution, the same exponent of the other two models 
(see 0). 

We present the results in Figures ^ to 0] in Figure ^ wc show E2 for 10 graphs 
on 100 nodes produced with BA model. In Figure[3we show E% for 10 7 graphs on 
100 nodes, produced with KR model, with r — 0.5. In Figure|3]we show E2 for 10 8 
graphs on 100 nodes, produced with DMS model. Oscillations in this case become 
more apparent using logarithmic scale, shown in Figure 

Here, the energy spectra has a maximum close to 100 3 , although the observed 
energy levels are far from this value in the examples. Moreover, there are no 
forbidden levels at least for E2 < 40, 000. This fact also shows that the oscillations 
are not caused by the discreteness nor the sparsity at high values of the energy. 
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BA, 100 nodes 
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Figure 1. Preferential attachment, 100 nodes, histogram over 10 7 graphs 

KR, 100 nodes, r=0.5 
16000 i , , , , 1 




120000 160000 200000 



Figure 2. Edge redirection, 100 nodes, r — 0.5, histogram over 
10 7 graphs 



The oscillations persist even varying the width of the bins in the histogram. 
This could be caused by the fractal pattern in the energy distribution, since each 
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DMS, 100 nodes 
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Figure 3. Attach to edges, 100 nodes, histogram over 10 s graphs 

oscillation exhibits a strong self-similarity with a sequence of sub-oscillations, see 
Figure |5] 



DMS, 100 nodes, log scale 
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Figure 4. Attach to edges, 100 nodes, histogram over 10 s graphs, 
logarithmic scale in Ei 
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KR, 100 nodes 
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Figure 5. Edge redirection, 100 nodes, r = 0.25, histogram over 
10 8 graphs 

It is interesting to note that the three models shown above share the same node- 
degree distribution. In all of the cases, degrees decay as a power-law of exponent 
3. However, the energy distribution presents some differences. For instance, it is 
clear that their modes differ: they are close to 4400 (BA), 5700 (KR), and 20000 
(DMS). Their oscillations also differ both in the place where they begin and in their 
amplitudes. 

Also, for r close to zero, the KR model gives high energy graphs, together with 
an interesting behavior of the oscillations, which become highly amplified. We con- 
jecture that a phase transition occurs for some specific value r = r c , a phenomenon 
similar to that in j^lj, which deserves further attention. We show in Figure the 
energy distribution for r = 0.25. 

The fractal pattern is clearly shown in the integrated density of energies, 



As was pointed out in [5] , the analysis of the cumulative distribution is better than 
binning in order to avoid statistical fluctuations. In FigureQ Figure|Sl and Figure 
13 we show different amplifications of N2(E) for KR with r — 0.25. Clearly, the 
structure resembles a Cantor devil staircase. 

2.5. Discrete Scale Invariance and log-periodic oscillations. Let us consider 
a model in which the first node has no outgoing edge, while the others have exactly 
one, as it happens in KR or BA models. Then, a graph with d nodes has d — 1 
edges, and for each edge x ~ y, \d{x) - d(y)\ < d - 2. Thus, E 2 (G) < E™ ax := 
(d— l)(d—2) 2 . This level of energy is reached with probability 1 in KR when r = 0. 
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KR, 200 nodes 
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Figure 6. Edge redirection, 200 nodes, r = 0.25, histogram over 
2 x 10 7 graphs 

KR, 200 nodes 
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FIGURE 7. Edge redirection, 200 nodes, r = 0.25, 10 7 graphs 

Now, consider the graph Gk depicted in figure^! It has 1 node at the center, k 
nodes with degree 2, k nodes "in a second row", attached to the previous k nodes, 
and d — 2k — 1 nodes with degree 1 and at distance 1 from the origin. 
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KR, 200 nodes 
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Figure 8. Edge redirection, 200 nodes, r = 0.25, 10 7 graphs 
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Figure 9. Edge redirection, 200 nodes, r = 0.25, 10 7 graphs 



The energy of Gk is given by 

e k : = E 2 {G k ) = (d - 2k - l)(d - k - 2) 2 + k(d - k - 3) 2 + k. 
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Figure 10. Graph with k edges in a second row 



If k < d, 

e k ~ (d-k- l) 3 . 

Within the KR model, if r is close to 0, this sort of graphs have non-negligible 
probability, and they explain the peaks at energies close to E™ ax , which can be 
seen in Figure ^2 If one pays a closer attention to one of these peaks, one will see 
that each of them has sub-peaks with smaller probability (see Figure 11211 . 
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Figure 11. Sub-peaks at high energies. KR, 200 nodes, r = 0.25, E 2 - 

These sub-peaks are due to graphs similar to but for which the k nodes at 
second row do not attach to those in first row so regularly. The number of energy 
levels around the peak associated to Gk is related to the number of graphs on k 
nodes (more precisely, to the number of forests, i.e. disjoint union of trees, on 
k leaves). This means that when k grows, these peaks get wider, and eventually 
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they become close enough to seem oscillations. Such a phenomenon can be seen in 
Figure El 
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Figure 12. Sub-peaks become close. KR, 2000 nodes, r = 0.20, E 2 . 
Since, for k <C d, 

ln(e fc )-Me fc+1 )~^f^~^, 

the peaks can be considered a log periodic structure. The log periodicity of oscil- 
lations was related to the existence of a discrete scale invariance (as opposed to a 
continuous one), which is due to the existence of preferred scaling ratios. For self 
similar fractals like the triadic Cantor set or the Sierpinski gasket, the existence 
of complex fractal dimensions were proposed, since the imaginary part leads to log 
periodic corrections of the scaling. 

In our setting, the existence of preferred scaling ratios correspond to the existence 
of different hierarchies in the distribution of node degrees in a given graph. That 
is, the degrees are not evenly distributed, but lacunary. 

We wish to stress that this fact was obtained independently for the degree distri- 
bution of hierarchical scale free graphs in |21| . and some models of graph generation 
were proposed. Here, the energy shows that those kind of networks are obtained 
also from classical models of graphs generation like the BA, KR, or DMS. 



3. Conclusions 

In this work we introduced the linking energy E 2 , and computed it for several 
graphs from 50 to 100,000 nodes. We generated the graphs with preferential at- 
tachment, edge redirection and attachment to edges models. These models are 
known as responsible of power law distribution of the nodes degree, and the graphs 
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generated are considered scale free. We also consider random (Erdos-Renyi) and 
regular lattice graphs. 

The frequency distribution of energies shows log periodic oscillations in the tail 
for the graphs generated with scale free models, which are associated to complex 
fractal dimensions and discrete scale invariance. They are not present for random 
or regular graphs. 

There are two interpretations of the terms scale free. In the first one, the range 
of the nodes degree distribution crosses several scales, instead of being concentrated 
around the mean number of connections. The second one suggest that the nodes 
degree distribution presents details at every interval of the range following a power 
law. 

Our results show the presence of discrete scale invariance in high energies graphs: 
they are scale free in the first sense, although the nodes degree are not evenly dis- 
tributed between the minimum and maximum values, but clustered in hierarchies. 
This is a fact observed in small networks gathered from real data. Also, for small 
energies, the distribution does not show the oscillations and the graphs seems to be 
scale free in the second sense. Hence, the models considered are able to generate 
both kind of graphs, and the energy could be used to discriminate between them. 
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